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ABSTRACT 

A refined, third-order plate theory that accounts for the 
transverse shear deformation is presented, the Navier solutions are 
derived, and its finite element models are developed. The theory does 
not require the shear correction factors of the first-order shear 
deformation theory because the transverse shear stresses are represented 
parabol ical ly in the present theory. A mixed finite element model that 
uses independent approximations of the displacements and moment 
resultants, and a displacement model that uses only displacements as 
degrees of freedom are developed. The mixed model uses C° elements for 
all variables and the displacement models use elements for the 
transverse deflection and C° elements for other displacements. 

Numerical results are presented to show the thickness effect on the 
deflections, and the accuracy of the present theory in predicting the 
transverse stresses. A comparison of the results obtained using the 
finite element models of the present theory with the experimental and 
the three-dimensional elasticity theory shcvs that the present theory is 
more accurate than the first-order shear deformation plate theory. 



1. BACKGROUND: REVIEW OF THE LITERATURE 

The classical laminate theory (CLT), which is an extension of the 
classical plate theory (CPT) to laminated plates, is inadequate for 
laminated plates made of advanced filamentary composite materials. This 
is because the effective elastic modulus to the effective shear modulus 
ratios are very large for such laminates. In addition, the classical 
plate theory is plagued with the inconsistency between the order of the 
differential equation and the number of boundary conditions (see Stoker 
[1]). An adequate description of the transverse shear stresses, 
especially near the edges, can be achieved with the use of a shear 
deformation theory. 

The shear deformation plate theories known in the literature can be 
grouped into two classes: (1) stress-based theories, and (2) 

displacement-based theories. The first stress-based shear deformation 
plate theory is due to Reissner [2-4]. The theory is based on a linear 
distribution of the inplane normal and shear stresses through the 
thickness, 

M 1 z M 2 z 

" (h 1 2 /6) < h/2 >’ 02 " (h 2 /6) < h / 2 >’ 

(1) 

where (c^.c^) and a g are the normal and shear stresses, (M-^, M 2 ) and Mg 
are the associated bending moments (which are functions of the inplane 
coordinates x and y), z is the thickness coordinate and h is the total 
thickness of the plate. The distribution of the transverse normal and 
shear stresses (a^, and o^) is determined from the equilibrium 
equations of the three-dimensional elasticity theory. The differential 
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equations and the boundary conditions of the theory are obtained using 
Castigliano 1 s theorem of Least Work. 

The origin of displacement-based theories is apparently attributed 
to Basset [5], who begins his analysis with the assumption that the 
displacement components can be expanded in series of powers of the 
thickness coordinate z. For example, the displacement component u^ 
along the x-direction in the N-th order theory is written in the form 

N , , 

Ui(x,y,z) = u(x,y) + z z n ^ n '(x,y) (2a) 

1 n=l x 

where x and y are the Cartesian coordinates in the middle plane of the 
plate, and the functions ii/ n ^ have the meaning 

d n u 

4 n) ( x ^)=— ^ , n = 0,1,2,... (2b) 

x dz n z=0 

Basset's work has not received as much attention as it deserves. In a 
1949 NACA technical note, Hildebrand, Reissner and Thomas [6] presented 
(also see Hencky [71) a first-order shear deformation theory for shells 
(which obviously can be specialized to flat plates). They assumed the 
following displacement field [a special case of Eq. (2a) for N = 1], 

u L (x,y,z) - u(x,y) + z *> x (x,y) 
u 2 (x,y,z) = v(x,y) + z i> y (x,y) (3) 

u 3 (x,y,z) = w(x,y) 

The differential equations of the theory are then derived using the 
principle of the minimum total potential energy. This gives five 
equilibrium equations in the five generalized displacement variables, 
u, v, w, to and q . 

a Jr 

The shear deformation theory based on the displacement field in Eq. 
(3) for plates is often referred to as the Mindlin plate theory. 
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Mindlin [8] presented a dynamic analysis of He.ncky's theory [7], and 
used the displacement field given by Eq. (3) for the vibration of 
isotropic plates. Historical evidence (from the review of the 
literature) points out that the basic idea of the displacement-based 
shear deformation theory came from Basset [5] and Hencky [7]. The shear 
deformation theory based on the displacement field given by Eq. (3) will 
be referred to as the first-order shear deformation theory. 

Following these works, many extensions and applications of the two 
classes of theories were reported in the literature (see [9-35]). 

Gol ' denveizer [16] has generalized Reissner's theory [2-4] by replacing 
the linear distribution of stresses through thickness [see Eq. (1)] by a 
distribution represented by an arbitrary function 4 >(z): 

= M 1 4>(z), o 2 = M 2 <}.(z), Og = Mg$(z) . (4) 

Kromm [14,15] presented a shear deformation theory that is a special 
case of Gol 'denveizer' s extension of Reissner's theory. In Kromm' s 
theory, the function <t>(z), instead of being arbitrary, is determined 
such that the transverse shear stresses vanish on the bounding planes of 
the plate. The displacement field in Kromm' s theory is of the form 


* H ♦ ! zo - 1 ^h x 


■I 


_ aw 3 7 n 4 z , 

U 2 ' U ~ Z W + 2 2 (1 - 3 7 )( V 


( 5 ) 


u 


3 


w 


where u, v, w, id and i|> are displacement functions which are functions 
x y 

of x and y only. Schmidt [29] presented an extension of Kromm 's theory 
by accounting for moderately large deflections (i.e., in the von Karman 
sense) . 
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Extension of the displacement-based theory to the moderately large 
deflections case is due to Medwadowski [17] and the extension to 
laminated plates is due to Whitney [23] and Whitney and Pagano [24]. 

The second- and higher-order displacement-based shear deformation 
theories have been investigated by Nelson and Lorch [27], Librescu [28] 
and Lo, Christensen and Wu [301. These higher-order theories are 
cumbersome and computationally demanding because with each additional 
power of the thickness coordinate an additional dependent unknown is 
introduced (per displacement) into the theory. 

Levinson [32] and Murthy [33] presented a third-order theory that 
assumes transverse inextensibility. The nine displacement functions 
were reduced to five by requiring that the transverse shear stresses 
vanish on the bounding planes of the plate. However, both authors (and 
also Schmidt [29]) used the equilibrium equations of the first-order 
theory in their analysis. As a consequence, the higher-order terms of 
the displacement field are accounted for only in the calculation of the 
strains but not in the governing differential equations or in the 
boundary conditions. Recently, Reddy [34,35] corrected these theories 
by deriving the governing differential equations by means of the virtual 
work principle. The theory presented in [34] accounts for moderately 
large rotations but is limited to orthotropic plates, while that in [35] 
is for the small-deflection theory of laminated plates. Finite element 
models of these theories were presented in [36,37] . The present report 
contains a summary of the research reported in references [34-37] on the 
development and analysis of a higher-order plate theory. 
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2. MATHEMATICAL FORMULATION OF A HIGHER-ORDER THEORY 


Consider a rectangular laminate with planform dimensions a and b 
and thickness h as shown in Fig. 1. The coordinate system is taken such 
that the x-y plane coincides with the mid-plane of the plate, and the z- 
axis is perpendicular to that plane (x^ = x, X£ = y and X 3 = z). The 
plate is composed of perfectly bonded orthotropic layers with the 
principal material axes of each layer oriented arbitrarily with respect 
to the plate axes. 

2.1 Displacement Field 

To obtain a parabolic distribution of the transverse shear stresses 
through the thickness, a cubic expansion in the thickness coordinate z 
[i.e., set N=3 in Eq. (2a)] is used for the inplane displacements and 
transverse inextensibility (i.e., the transverse normal strain is 
zero) is assumed. The additional displacement functions are determined 
in terms of the five displacement functions in Eq. (3) by requiring the 

u. 

transverse shear stresses to vanish at z = ± 2 • The resulting 
displacement field is given by 

u^x.y.z) = u(x,y) + z[* x (x,y) - | (|] 2 (^ x (x,y) + (x,y))] 

u 2 (x,y,z) = v(x,y) + z[* y (x,y) - | (f] 2 O y (x,y) + fy (x,y))] 

u 3 (x,y,z) = w(x,y) , (6) 

where u^, U 2 and U 3 are the displacements in the x-, y- and z- 

directions, respectively. The displacements of a point (x,y,0) on the 

midplane of the plate are denoted by u, v and w; ip and are the 

x y 

rotations of the normals to the midplane about the y- and x-axes, 
respectively. The assumed deformations of the transverse normals in 



various plate bending (displacement based) theories are shown in Fig. 

2 . 

The strains of the von Karman theory (i.e., only the products and 
squares of aw/ax and aw/ay in the strain-displacement equations of the 
large deflection theory are retained) can be obtained from Eq. (6): 



2.2 Stress Field 


Assuming that each layer of the laminate possesses a plane of 
material symmetry parallel to the x-y plane, the constitutive equations 
for the k-th layer can be written as 
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where a. and e. (i = 1,2, 4, 5, 6) are the stress and strain components 

referred to lamina coordinates and Q. .'s are the plane-stress reduced 

^ J 

elastic constants in the material axes of the k-th lamina. 


Qn E 1 /(1 -v 12 m 21 ) . Q 22 Q u , Q 12 ^ 12^22 

G 44 = G 23 * ^55 = G 13 ’ G 66 = G 12 


and 






and 


G. . 
TJ 


are the usual engineering constants. 


2.3 The Principle of Virtual Displacements: Governing Equations 

The equations of equilibrium can be obtained using the principle of 
virtual displacements. In analytical form, the principle can be stated 
as follows (see [38]): 

J (fill + <sW)dV + SV = 0 (11) 

0 


where U is the total strain energy density due to deformation, W is the 

potential of distributed external loads, V is the potential of discrete 

external loads, and n is the volume of the laminate and 6 denotes the 

variational symbol. Integrating the expressions in Eq. (11) through the 

laminate thickness, integrating by parts with respect to x and y, and 

collecting the coefficients of <su, 6v, 6w, 5 ^ , and 6ii> , the following 

x y 

equations of equilibrium are obtained: 


su: N. + N, = 0 

— l,x 6,y 


sv: N, + N„ = 0 

— 6,x 2,y 


Sw: Q, +Q„ - — — (R, +R 0 1+ — — — f P . +2P c +P 0 ) = q 

— H l,x x 2,y h2 1 1 , x 2,y J 7h2 k l,xx 6,xy 2,yy J 


3h 2 
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4 


6 V 

Mi +M C - 
1.x 6,y 

Q l + 

4 

h 2 

Sii> : 

M c +M„ - 

q„+ 

4 


6,x 2,y 

x 2 

h 2 


h 2 x 3h 2 

_ _4_ 
h 2 3h 2 


(P. +P C 1=0 
v 1.x 6,y ; 


( 12 ) 


where 


h/2 3 

(N . ,M . ,P . ) = J a.(l ,z,z^)dz (1=1,2, 6) 

1 1 1 -h/2 1 

h/2 ? 

(Qo.R ? ) = I ° 4 (1,z )dz 

h/2 ? 

(Qi.Ri) = J ° 5 0,z )dz 

1 1 -h/2 5 


(13) 


2.4 Laminate Constitutive Equations 

The laminate constitutive equations are relations between the 
stress resultants (13) and the strains (7). Using Eqs. (7) in (9) and 
the result in Eq. (13), we obtain 
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where A^j, B ^ , etc., are the laminate stiffnesses, defined by 


h/2 p o n c 

(» ii ,fl ii ,O ii ,E n ,F,,,H ji )=J Q (l,z,zSz\z\z D )dz (1,j 
>j ij iJ ' J 'J 1 J -h/2 4 


(A. .,0. . ,F . .) 
v ij* ij -\y 


h/2 P 4 

J Q , .(l,z ,z 4 )dz (i , j = 4,5) 
-h/2 


(15) 

= 1 , 2 , 6 ) 

(16) 
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3. EXACT CLOSED-FORM SOLUTIONS 


In the general case of arbitrary geometry, boundary conditions and 
lamination scheme, exact analytical solutions to the set of differential 
equations in Eq. (12) cannot be found. However, closed-form solutions 
for the static case exist for certain ‘simply supported* rectangular 
plates with two sets of laminate stiffnesses (see Reddy [34,35]), as 
described below. 

It is possible to obtain the Navier-type solutions for the 
following two types of simple supported boundary conditions. 

u(x,0) = u(x,b) = v(0,y) = v(a,y) = 0 
S-l : (17a) 

N 2 (x, 0) = N 2 (x,b) = N^(0,y) = N^a.y) = 0 

u(0,y) = u(a,y) = v(x,0) = v(x,b) = 0 
S-2: (17b) 

N 6 (0,y) = N 6 (a,y) = N & (x,0) = Ng(x,b) = 0 

w(x,0) = w(x,b) = w(0,y) = w(a,y) = 0 

S-l P ? (x,0) = P ? (x,b) = P-, (0,y) = P,(a,y) = 0 

and : ^ 1 1 1 (17c) 

S-2 M 2 (x,0) = M 2 (x,b) = M^O.y) = M^a.y) = 0 

4 x (x, 0 ) = ip x (x,b) = * y (0,y) = iP y (a,y) = 0 
Here a and b denote the planform dimensions along the x and y 
coordinates, respectively, and the origin of the coordinate system is 
taken at the lower left corner of the rectangular plate (see Fig. 1). 

The Navier procedure involves expressing u, v, w, , and in 

* y 

terms of the Fourier series with undetermined coefficients. The 
functions in the series are selected such that they satisfy the boundary 
conditions of the plate. The coefficients are then determined by 
satisfying the equilibrium equation, Eq. (12). 
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Exact solutions can be obtained for rectangular plates with the 
combination of boundary conditions given in Eq. (17) and with the 
following zero plate stiffnesses [see Eqs. (14, (15) and (16)]: 

For S-l boundary conditions [Eqs. (17a) and (17c)]: 


A 

E 

A 


16 

16 

45 


" A 26 " B 16 
= E 26 = F 16 


°45 F 45 


= B 26 ' °16 


F 26 H 16 


= D 


26 


= H 


26 


= 0 
= 0 


(18a) 


For S-2 boundary conditions [Eqs. (17b) and (17c)]: 


A 16 A 26 B 11 B 12 °16 °26 0 
E 11 = E 12 = F 16 = F 26 = H 16 = H 26 = ° 

(18b) 

A 45 = °45 = F 45 = 0 



The stiffnesses of the general cross-ply laminate satisfy Eq. (18a) 
while the stiffnesses of the antisymmetric angle-ply laminate satisfy 
Eq. (18b). A complete discussion of the exact solutions is presented in 
[38]. 
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4. FINITE-ELEMENT MODELS 

4.1 Displacement Model 

Although the present higher-order theory has the same five 

generalized displacements as the first-order theory, second-order 

derivatives of the transverse deflection w appear in the total potential 

energy expression of the higher-order theory. An examination of the 

total potential energy functional of the higher-order theory reveals 

that u, v, w, ip , i|) , and and appear in the geometric boundary 

conditions. Therefore, the finite element model based on the total 

potential energy requires the C°-continuity of u, v, i> and i> and the 

x y 

(^-continuity of w across interelement boundaries. Thus, the Lagrange 

interpolation of u, v, and ip , and the Hermite interpolation of w are 

x y 

required. Consequently, the classical plate bending (i.e., Hermite) 

element can be obtained as a special case from the present element by 

suppressing i|i and i|> degrees of freedom. In the present study, the 
x y 

four-node rectangular element with u, v, w, ip , i|i , and as the 

x y 3x ay 

nodal degrees of freedom is used. 

The elemental equations of the displacement model are of the form 
(see [36]) 

[K e ]{4 e } - |F e ) (19) 

where {A e [ denotes the vector of generalized nodal displacements 
and (F e | is the force vector that contains the applied transverse load 
as well as the contributions from the boundary of the element (see Reddy 
[39]). 

4.2 Mixed Model 

To relax the continuity requirements placed in the displacement 
model described above (i.e., to reduce the C^-continuity to the C°- 
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continuity on w), a mixed formulation of the problem is considered. In 
a mixed model, independent approximations of the displacements and 
stress resultants are used. A close look at Eq. (12) shows that a mixed 
variational formulation of these equations can be obtained by treating 
u, v, w, i|> x , ip , M^, M 2 , Mg, P^, P 2 and Pg as the nodal degrees of 
freedom. The governing equations for u, v, w, i|» x and ii» are given in 
Eq. (12). The equations for the other six variables (M-j, M 2 , Mg, P-j, P 2 
and Pg) are provided by the laminate constitute equations (14) and (15). 

The mixed model is of the form (see [37]) 

[K e ]{A e ] = {F e } (20) 

where [K e ], {a 6 } and {F e } are the generalized element stiffness matrix, 
element displacement vector and element force vector, respectively. 
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5. NUMERICAL RESULTS 


Numerical results obtained using the displacement and mixed finite 
element models described herein are presented for the bending of 
laminated anisotropic plates. Only respresentative results are included 
here. For additional results, see Refs. 34 thru 37. 

In all the problems considered, the individual layers are taken to 
be of equal thickness, and only two sets of material properties are used 
(these values do not necessarily indicate the actual properties of any 
material, but they serve to perform parametric studies): 

Material 1 : 

E 1 /E 2 =25, G 12 /E 2 =0.5, G 23 /E 2 =0.2, v 12 =0.25 (21) 

Material 2 : 

E^E^O, G 12 /E 2 =0.6, G 23 /E 2 =0.5, v 12 =0.25 (22) 

It is assumed that G^ = G-^, v^ 3 = v^. For the analysis based on the 
first-order shear deformation plate theory (FSPT), unless stated 
otherwise, the shear correction coefficients are taken to be 

= K 2 2 = 5/6. 

First, the (mixed) element was evaluated for its sensitivity to 
locking by using different integration rules. A two-layer, cross-ply 
[0/90], simply supported square plate under uniformly distributed 
transverse load was analyzed using the uniform 4x4 mesh of linear 
elements and 2x2 mesh of quadratic elements (in a quadrant), and various 
types of integration rules. Full integration refers to the usual 
integration rule (2x2 for the linear element and 3x3 for the quadratic 
element) for evaluating both bending and transverse shear terms of the 
stiffness matrix. Reduced integration refers to one point less (in each 
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coordinate direction) than that used in the usual integration scheme. 
Mixed (selective) integration refers to the use of the full integration 
rule for bending terms and the reduced integration rule for the shear 
terms. From the numerical results presented in Table 1, it is clear 
that full integration for very thick plates (4 < ^ < 10), mixed 
integration for moderately thick plates (10 < ^ < 25) and reduced 
integration for thin plates (50 < ) give the best results. Also, the 

quadratic element is less sensitive to integration rules than the linear 
element. In general, mixed integration gives the best results for the 
entire range of thicknesses. 


5.1 Thickness Effect 

The accuracy of the present higher-order theory over the first- 
order and classical plate theories is demonstrated using a specific 
problem. Consider a four-layer [0/90/90/0] square cross-ply laminated 
plate (Material 1) subjected to sinusoidally distributed transverse 
loading. The plate is simply supported along all four edges with the 
in-plane displacements, unconstrained in the normal direction, and 
constrained in the tangential direction. The deflections and stresses 
are non-dimensional i zed as follows: 


a a 100E 2 h ' 
w = w^, °) 4~~ 

V 


— ,a a hv _h_ 

°i - °]^2» 2 s 2 


V 


°2 

°6 


, ^2 
/a a hv h 

a 2''2’ 2’ 4^ 2* 

V 

° 6 (0, °* ~~ 2 ’ 

V 


(23) 
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a 5 = o 5 (0, f, 0) h 


q a* 
H o 


= a 4 (§, 0, 0) 


V 


where q Q denotes the intensity of the distributed transverse load, 
q(x,y) = q Q sin sin The exact maximum deflections obtained by 
various theories are compared for various side-to-thickness ratios in 
Fig. 3. The deflections predicted by the present higher-order shear 
deformation theory (HSOT) are in excellent agreement with the 3-D 
elasticity solutions of Pagano and Hatfield [40] . In the range of 
medium to thick plates (5 < ^ < 20), the HSDT yields more accurate 
results than the first-order shear deformation theory (FSDT). The 
effect of shear deformation on the deflections is apparent from the 
difference between the deflections predicted by the classical laminate 
theory (CLT) and the two shear deformation theories. The shear 
deformation causes a reduction in stiffness of composite plates and this 
reduction increases with decreasing side-to-thickness ratio. Clearly, 
the CLT is adequate in predicting the deflections of thin laminates 
(a/h > 50). 


5.2 Throuqh-the-Thickness Distribution of Transverse Shear Stresses 
The exact stresses a 4 e and e a^ z computed using the 
constitutive equations of the HSDT are greatly improved over the results 
obtained using the CLT (classical laminate theory) and FSDT (first order 
shear deformation theory) as shown in Fig. 4a. The shear stresses 
obtained using the constitutive equations are on the low side of the 3-D 
elasticity solutions. This error might be due to the fact that the 
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stress continuity across each layer interface is not imposed in the 
present theory. As in the case of the classical laminate theory (CLT), 
the transverse shear stresses can also be determined by integrating 
equilibrium equations (of three-dimensional elasticity in the absence of 
body forces) with respect to the thickness coordinate: 


xz 


h/2 

I (°v v + °v» J dz 

-h/2 


x,x xy,y' 


h/2 

a - - f (a +o )dZ 
yz J _ h/2 V xy,x y,y ; 


(24) 


The foregoing approach not only gives single-valued shear stresses at 
the interfaces but yields excellent results for all theories in 
comparison with the three-dimensional solutions. Despite its apparent 
advantage, the use of stress equilibrium conditions in the analysis of 
laminated plates is quite cumbersome. Typical stress distributions of 
erg = a xz and o q = a through the thickness (a/h = 10) are shown in 
Fig. 4b. Note that the stress discontinuity at the laminate interfaces 
is due to the mismatch of the transformed material properties. 


5.3 Relative Magnitudes of Transverse Stresses 

The orders of magnitude of the stresses can be assessed by 
comparing the solutions of the higher-order theory with those of the 
three-dimensional theory of elasticity. To this end, an analytical 
solution is developed using the plate equations and boundary conditions 
in three dimensions following the procedure of Pagano [40]. The 
results, for specific problems, are shown in Tables 2 and 3. 

The ratios of stresses in a simply supported square isotropic plate 
under sinusoidally (SDL) and uniformly (UDL) distributed transverse 
loads are given in Table 2. The number of terms used in the double- 
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TABLE 1 Effect of Reduced Integration and Type of Mesh on the Accuracy 
of the Solution (Cross-Ply [0/90] Simply Supported Square 
Plate Linder UDL; Material 1]) 


Type of 


Type of Mesh** 

<f> 

Type of Mesh 

Integration 

4x4L 

2x2Q 

4x4L 

2x2Q 

Full 


3.0792 

6.9705 


1.4784 

1.7479 

Mixed 


3.1045 

3.0723 

20 

1.7475 

1.7534 

Reduced Sol. 

4 

3.1145 

3.0769 

1.7576 

1.7548 

Exact 


3.0706 

3.0706 


1.7509 

1.7509 

Full 


2.5728 

2.5813 


1.3261 

1.7252 

Mixed 


2.5990 

2.5808 

25 

1.7270 

1.7334 

Reduced 

5 

2.6090 

2.5850 

1.7371 

1.7348 

Exact Sol . 


2.5791 

2.5791 


1.7310 

1.7310 

Full 


1.8409 

1.9192 


0.2826 

1.6789 

Mixed 

10 

1.9185 

1.9193 

100 

1.6927 

1.7001 

Reduced 

1.9288 

1.9218 

1.7029 

1.7014 

Exact Sol . 


1.9173 

1.9173 


1.6977 

1.6977 


w = wE^h^ 100/q^. 

4x4L = 4x4 uniform mesh of linear elements 
2x2Q = 2x2 uniform mesh of quadratic elements 




TABLE 2 Stress Ratios in a Simply Supported Square Isotropic Plate 

Under Different Loadings (SSL-Sinusoidal Loading, UDL-Uniformly 
Distributed Loading) (E = 10° psi, v = 0.3, q Q = 1.0 psi.) 


(-) 

V 

SSL 

UDL (m = 

n = 30) 

O/O 

Z X 

a / 0 

xz' xy 

0 / 0 
T X 

a /a 
xz' xy 

4 

0.28740 

0.53731 



5 

0.19069 

0.43579 

0.13338 

0.47597 

6.25 

0.12478 

0.35207 

0.08667 

0.38747 

10 

0.04989 

0.22262 

0.03443 

0.24861 

12.5 

0.03210 

0.17859 

0.02212 

0.20048 

20 

0.01261 

0.11197 

0.00868 

0.12654 

25 

0.00808 

0.08964 

0.00556 

0.10148 

50 

0.00202 

0.04486 

0.00139 

0.05092 

100 

0.00051 

0.02244 

0.00035 

0.02548 


Note: o = a and a = a 

x y xz yz 

a ,o ,a values are at (x = y = -|, z = ■!}) 
x y z v J 2 2' 

a xz is at (x = 0, y = z = 0) 

a yz is at (x = ^ y = 0, z = 0) 

L. 

a xy is at (x = y = 0, z = ^ 


TABLE 3 Stresses in a 4-ply [0/90/90/0] Square Plate Under 
Sinusoidal Loading: Material 1 (q Q = 1.0 psi.) 


(£> 


°x 

-a 

y 

O 

z 


a 

XZ 

1 

Q 

<< 

N 

°xy 

4 


11.524 

45.240 

1.0 
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Fourier series, for the UDL case, is m=n=30, and for SDL, we take m = n 

= 1. It is clear that the ratio of normal stress a z to the inplane 

normal stresses a and a is very small (less than 5%) for moderately 
x y 

thick plates and very thin plates i.e., (10 < ^ < 100). However, the 

transverse shear stresses o* z and o yz are about 25% of the inplane shear 

stress, a , for moderately thick plates (10 < r- < 25). In Table 3, 
xy n 

similar results for a 4-ply [0/90/90/0] laminate are presented (Material 
1 and sinusoidal loading); the same trend as described for isotropic 
plates is seen, confirming the fact that the transverse shear stresses 
should be considered before considering the transverse normal stress for 
moderately thick and thin plates. 

5.4 Geometric Nonlinearity Effect 

Next, the results of the von Karman (nonlinear) theory are 
presented. In Figures 5 and 6 , the mixed finite element results of the 
higher-order theory are compared with the experimental and thin plate 
theory results of Zaghloul and Kennedy [41,42]. In Figure 5, center 
deflection versus the load intensity are presented for a simply 
supported, orthotropic, square plate (a = b = 12 in. , h = 0.138 in.) 
under uniform loading. The material properties used are E-| = 3x10® psi, 
E£ = 1.28x10® psi, G ^2 = = G 23 = 0.37x10® psi, v = 0.32. The finite 

element results are in excellent agreement with the experimental 
results. The thin plate theory results [41] are in considerable error, 
which can be attributed to the neglect of transverse shear 
deformations. In Figure 6 , the center deflection of a clamped 4-ply 
symmetric bidirectional [0/90/90/0] square plate (a = b = 12 in, h = 
0.096 in) under uniform loading is presented. The material properties 
used are: E-| = 1.8282x10® psi, E 2 = 1.8315x10® psi, G^ = G^ = G 23 = 



.3125xl0 6 psi, \j = 0.23949. The finite element results obtained using 
the higher-order theory are almost the same as those obtained using the 
first-order theory, and are in better agreement with the experimental 
results than those obtained using the classical theory. The difference 
between the experimental and higher-order theory solutions is attributed 
to inexact simulation of the experimental boundary conditions (clamped), 
which have significant effect on the deflections in composite laminates 
(because of the degree of orthotropy) than in homogeneous orthotropic 
plates. 
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6. SUMMARY AND CONCLUSIONS 

A higher-(third-)order shear deformation theory is presented and 
its accuracy is demonstrated by comparing with the 3D elasticity theory, 
the classical plate theory and the first-order shear deformation 
theory. Mixed and displacement finite element models of the theory are 
described. The finite element solutions for nonlinear bending are 
compared with experimental solutions and found to be in good agreement. 

The higher-order shear deformation theory gives, in general, more 
accurate solutions for the bending of laminated anisotropic plates than the 
classical plate theory and the first-order shear deformation theory. The 
present observations concerning the effects of shear deformation, and 
coupling and material anisotropy on the bending of plates are in conformity 
with other investigators' findings (see [23,24]). The first-order shear 
deformation theory (FSDT) is found to overpredict the deflections, and to 
underestimate the natural frequencies and buckling loads (results not 
included here) for medium to thick anti -symmetric cross-ply and angle-ply 
plates. In FSDT, the shear correction factors play a crucial role in the 
determination of deflections. No such corrections are required in the 
present higher-order theory. In addition, the finite element based either 
on the displacement formulation or the mixed formulation is less 
susceptible to locking than the element based on the displacement 
formulation of the first-order theory. 

The present theory gives accurate through-the-thickness distribution 
of the interlaminar shear stresses but does not contain the transverse 
normal stress. It is possible to account for a by expanding w as w(x,y) = 
w Q (x,y) + zii» 2 (x,y), which would increase the number of nodal degrees of 
freedom by one. The discontinuity of the stresses (obtained using 
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constitutive equations) at layer interfaces can be eliminated by imposing 
continuity requirements, which might lead to a more complicated 
computational model. 
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Kinematics of various plate theories 





UNDEFORMED 


DEFORMED IN CLASSICAL 
(kirchhoff) THEORY 


DEFORMED IN THE FIRST 
ORDER THEORY 


DEFORMED IN THE THIRD 
ORDER (PRESENT) THEORY 


Figure 2. Assumed deformation patterns of the transverse normals 
in various displacement-based theories. 
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(a) From constitutive equations (b) From equilibrium equations 


Fiqure 4 The distribution of normalized transverse shear stresses 
through the thickness of a four-layer [0/90/90/0] square 
laminated plate under sinusoidal load (a/h = 10). 
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CENTER DEFLECTION (in,) 



V \ I I l 1 I 1 1 1 1 

0 0,4 0,8 1,2 1.6 2,0 


LOAD INTENSITY (psi) 


Figure 5. Center deflection versus load intensity for a 
simply supported square orthotropic plate 
under uniformly distributed transverse load. 

The following simply supported boundary condi- 
tions were used: 

v = w = \jj - M = P = 0 on side x = a 
y x x 

u = w = = My = = 0 on side y = to 
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CENTER DEFLECTION (in,) 



it t i i i > i t i » i 

0 0,4 0,8 1,2 1,6 2,0 


INTENSITY OF TRANSVERSE LOAD (psi) 


Figure 6. Center deflection versus load intensity for a 

clamped (CC-1) square laminate [0 /90 /90 /0 ] 
under uniform transverse load (von Karman theory). 

CC-1: u = v = w = ip = i|> - 0 on all 

* y 

four clamped edges 
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